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Abstract 

A general approach for constructing multidimensional quasi-exactly 
solvable (QES) potentials with explicitly known eigenfunctions for two 
energy levels is proposed. Examples of new QES potentials are pre- 
sented. 
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1 Introduction 

The importance of exactly solvable potentials in quantum mechanics is well 
known. But even in one-dimensional case the number of such potentials 
is limited. Therefore, much attention has been given to the quasi-exactly 
solvable (QES) pontentials for which a finite number of energy levels and 
corresponding wavefunctions are known in explicit form. The first examples 
of such potentials were given in [TJ EJ HI E| and subsequently many QES 
potentials were established [SJ El M UH1 EEU EH EH] (for review see book by 
Ushveridze |Hj). 

In multidimensional case many exactly solvable systems admit separation 
of variables and thus they are reduced to one-dimensional problems [To] . 
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The investigation of multidimensional and many-body systems which are 
not amenable to separation of variables is important from mathematical and 
physical points of view. 

One of the most remarkable exactly solvable models in iV-body quan- 
tum mechanics is the Calogero model [TIH ITTj . This model describes one- 
dimensional iV-body problem with quadratic and inverse square interacting 
potential. Soon after these papers Sutherland extended Calogero model to 
a model where interaction takes place on a circle jTSl EH 12*0] . The further 
progress was connected with Olshanetsky-Perelomov integrable systems [21] 
(see also review [22] )• In [22] explicit examples of QES iV-body problems on 
a line were introduced for the first time. Later QES iV-body problems were 
studied in [21123 • 

Recently, there has been achieved some progress in studying exact solv- 
ability of such systems. In [20] a new class of QES many-body Hamiltoni- 
ans in arbitrary dimension was constructed. The multidimensional Darboux 
transformation was proposed in [2Z] and new examples of QES multidimen- 
sional matrix Schrodinger operators were presented there. The authors of 
[2H] have developed a systematic procedure for constructing exactly solvable 
and QES many-body potentials by purely algebraic means. The method of 
multidimensional supersymmetric (SUSY) quantum mechanics was applied 
to the investigation of iV-particle systems and an explicit construction of ex- 
actly solvable 3-particle as well as QES iV-particle problems on a line were 
presented [2H1 HUB E2[ • Very recently, two new methods based on higher-order 
SUSY quantum mechanics for the investigation of two-dimensional quantum 
systems, whose Hamiltonians are not amenable to separation of variables, 
were presented in [32] ■ 

In this paper we develop a simple method for generation of QES potentials 
with two known eigenstates in arbitrary dimension. For one-dimensional case 
the problem of constructing QES potential with two known levels was solved 
completely in the frame of SUSY quantum mechanics [33J EH E2] or using a 
simple method in which two wave functions are chosen in such a way that 
they lead to the same potential [3EIIHZ!- Just this method will be generalized 
for multidimensional case. 

Note that in fact the idea of the papers [33] EU HUH HUE EH! * s Dase d on the 
inverse method. For the first time this method was used by Ushveridze fHJ 
for construction QES potentials. It rests on the very simple idea: instead of 
looking for solution of Schrodinger equation for a given potential one should 
reconstruct the potential starting with appropriately chosen eigenfunctions. 
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We would like to mention also paper [3H] where detailed comparison of su- 
persymmetric approach proposed in and Turbiner-Ushveridze approach 
[H] based on finite-dimensional representations of sl(2,R) was done. 



2 Construction of QES multidimensional po- 
tentials 

We consider the Schrodinger equation 

Hip = Ei/j (1) 
in n-dimensional case with the Hamiltonian 

H=-^A + V{xt,---x n ), (2) 

where A = (V, V) = J^zj V is the nabla operator in n-dimensional space. 

The wave function of the ground state with the energy Eq has no zeros 
and can be written in the following form 

V>o = e~ F , (3) 



- n ■ 



where F is a nonsingular function of x±, . . . x Ti 

Substituting (J5J) into ((TJ) we express the potential in terms of the wave 
function ip : 

V = E + = E + \ \{VFf - AF] . (4) 

Note that this transformation of Schrodinger equation from linear form 
to non-linear one is well known in literature. See, for instance, review by 
Turbiner, where this procedure was used for construction of convergent per- 
turbation theory in quantum mechanics |39j . 

The wave function of the excited state with the energy E\ we write as 
follows 

fa = 0e" F , (5) 

where = ipi/ip®. This eigenfunction must lead to the same potential V and 
two eigenf unctions satisfy the following equation 

Eo + l^l = El + l^l. (6) 
2 A 2 Vi 
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In terms of and F functions this equation reads 

2(VF, V0) = A0 + 2e0, (7) 

where 

£ = Ei — E . 

In one-dimensional case this equation can be easily solved with respect to 
F |36| 13 7j. Then choosing different generating functions <fi we obtain different 
QES one-dimensional potentials with two known eigenstates. 

In multidimensional case the solution of equation (J7J) is a nontrivial task. 
In this paper we find the solution of (J7J) with respect to F for some special 
cases of function <p. Note that equation (JJJ) is a nonuniform linear equation 
with respect to F. Any solution of it can be written in the following form 

F = f + f, (8) 

where / is a particular solution of (J7|) and / is a general solution of the 
uniform equation 

(V/,V0) = O. (9) 

Direct substitution of function (JSJ) into equation (J7|) shows that F is indeed 
a solution of the equation. In order to satisfy square integrability of wave- 
functions (j3J) and (J5J) we have to choose such / and / that F — ► +oo when 
|x| — > +oo. 

Note that in one-dimensional case this equation gives only trivial solu- 
tion / = const, which does not influence the potential (j3J). New point of 
multidimensional case is that uniform equation has many nontrivial so- 
lutions /. This gives more possibilities for constructing QES potentials in 
multidimensional case for given in comparison with the one-dimensional 
case. 

Case 1. This case corresponds to the following form of <fi 

= f>(zi). (10) 

i=i 

Then equation (J7j) can be rewritten as 

n or n 

2 E^V = E(^ ,/ + 2 ^), (ii) 

i=l 1 i=l 
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where # = fj and cj>» = 



The particular solution of equation (JTTJl can be found in the form of the 
function with separated variables 



i=i 

where 



/ = £/<(**), (12) 

i=i 



dxi 2<fi' i 

and the constants Ai, . . . A n have to satisfy the following condition 



The solution of the uniform equation 

dx 



(13) 



E^ = o- (14) 



t<f>'S = (15) 



i=i 
reads 

/ = f(,XiM - X2(x 2 ), XiM - Xafa), Xi(xi) ~ XnM), (16) 
where / is an arbitrary function of n — 1 arguments and 

*M-Fm (17) 

Finally, the solution of (|TT|) reads 
F = l^ ^-771 + /(Xi(^i) - W2), • • •)• (18) 

Note that in the case of / = the variables of potential (JH) are separated. 
It is the function / that is responsible for nonseparability of variables in 
potential V. 

Case 2. Let us choose the generating function in the following form 



rw*i)- (1 9 ) 
1 
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i=i 



Now equation J7J) reads 



2E 



n 




(20) 



i=i 



<t>i dxi 



The solution is similar to the solution of the previous case. We can write it 



as 



i=l 



n 



I 



Xi + (~f + AQ0j(z) 
2$(z) 



d» + /(Xl(^l) - X2(^2), • • •) 



(21) 



where the constants Ai, A 2 , . . . also satisfy condition (JHJ), / is an arbitrary- 
function of n — 1 arguments and 



3 Examples 

QES potential is given by equation (@|, where the function F is represented 
by expression (J18)) in the case 1 and by expression (J21j) in the case 2. For 
this QES potential we know two energy levels E and E\ as well as the 
corresponding wave functions (0), ©• Choosing various functions (pi for the 
case 1 or the case 2 we obtain different QES potentials. Without loss of 
generality we choose in all expressions Eq = and then E\= e. 

The functions <pi and the parameters Aj must be chosen in such a way 
that the function / is a nonsingular one and ipo, ip\ are square integrable 
functions. Note that in multidimensional case we have also possibility to 
choose different / for this purpose. 

Example 1. Let us choose 



which corresponds to the case 1. For this generating function we obtain 





= o E a i x 1 



1 ^~ 



i=l 




a-i + Aj + diExl 
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To satisfy the nonsingularity of the functions /* we have to choose A, = —a^. 
Note that a, satisfy the same condition as Aj, namely = 0. Then 

fi{%i) T*^i ' lllXj. 

4 a. 



And finally, we obtain 



F(x x ,x 2 ,...) = ~;J2 x2 i+ f 



l/ai 1/ai 
£ 2 i I I x l X l 



4 4rf 2 M -V^ ' l/a 3 • 

here we rewrite f(— lnxi — - lnxo, — Inx, — - \nx-i, . . .) as some new func- 
tion 



/ 



1/aa ' 1/oa ' " " " 



-2 x 3 

Let us apply these results to the two-dimensional case. Now, a\ + a 2 = 
and we may choose a x = —a 2 = 1. Then 

QES potential and wave functions with zero and e energy levels read 



V(x,y) = - 



Cff - f" + j 



[x 



y 2 ) + exyf - - 



fa{x,y) = Cl (x 2 - 2/ 2 ) e -f(- 2 +^ 2 )-/>f), 



where Cq and C\ are normalization constants. 

Choosing / as some polynomial we reproduce the two-dimensional po- 
tentials studied in [10] as interesting examples of bottomless potentials with 
bound states. 

For / = we obtain the isotropic harmonic oscillator. Then, in Dirac 
notation, ip = |0, 0) is the ground state eigenf unction, and ipi — ^75 (|2, 0) — 
1 0, 2) ) corresponds to the second excited energy level. 

Note that ^0 corresponds to the ground state if / is a nonsingular func- 
tion. For singular / wave function ipo may correspond to an excited state. 
For example, let us choose / = —\n(xy). Then 

il) Q (x,y) = c xye~^ {x2+y2) 
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corresponds to the second excited energy level. 

Example 2. This example illustrates the case 2. We choose 



(p = Y[xi. 



i=l 



Using (|2*T|) we obtain 



^ = Ef|e + A i )^ + /(x?-4^-4-); 



i=i 



remember that X) A, = 0. 

Let us consider a particular three-dimensional case with Ai = A2 = A3 = 
and / = a(2x 2 — y 2 — z 2 ) 2 . Then 

v = (S " i6a ) x2 + " 4a ) {y2 + z2)+ 

+Aa{2x 2 -y 2 - z 2 ) 2 [2a{Ax 2 + y 2 + + |) - ~e, 



^0 = Coe -^^-^ 2 -^- 2 ) 2 5 



where Co and ci are normalization constants. 

For a = we have harmonic oscillator. Then t/> = |0, 0, 0) is the ground 
state function, and ipi = 1 1,1,1) corresponds to the third excited energy 
level. For a 7^ we obtain nontrivial QES three-dimensional potential, the 
variables of which cannot be separated. 

Here we would like to underline that question about separation of vari- 
ables in multidimensional Schrodinger equation is an interesting and non- 
trivial task (see, for instance, [MJII2])- We plan to discuss this problem in a 
part concerning our approach in a separate paper. 

4 Conclusions 

We developed a simple approach for constructing QES multidimensional po- 
tentials with two known energy levels and corresponding wave functions. The 
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proposed method is a direct extension of the general approach proposed ear- 
lier for one- dimensional case The central point of our approach is 
equation (j7J and the main problem is to solve this equation with respect to 
the function F. In contrast to one-dimensional case, when the corresponding 
equation can be easily solved, there is a nontrivial task to solve it in mul- 
tidimensional case. In this paper we find general solutions of equation (J7J 
and construct new multidimensional QES potentials for some special cases 
of function 0. Finding of other solutions of equation J7J is an interesting 
problem for further investigations. 

Let us stress the following new point which appears in multidimensional 
case in contrast to one- dimensional case. Namely, in one- dimensional case 
/ = const and therefore the ratio of the two eigenfunctions ipi/^o = 
and distance between corresponding energy levels e entirely determine the 
potential energy V j34[ 13^1 137] . In multidimensional case the potential energy 
V is determined by <f), f and e. Equation Q allows many nontrivial solutions 
for /. As result we obtain a family of potentials V for fixed <fi and fixed energy 
gap e choosing different solution for /. This feature of multidimensional case 
is explicitly shown in the example 1, where potential V(x, y) for fixed 
and e depends on arbitrary function f(xy). Thus, in multidimensional case 
there are more possibilities for constructing QES potentials for a given in 
comparison with the one- dimensional case. 

In one-dimensional case all eigenfunctions can be easily ordered using 
oscillation theorem which states that the n-th eigenfunction has n zeros. 
For multidimensional case it is known that ground state eigenfunction also 
has no zero and the corresponding energy level is non-degenerated [13], but 
there is no theorem which connects zeros and order of excited eigenfunctions. 
Therefore, one cannot give an exact answer between which states the energy 
gap is calculated. We can only state that if an eigenfunction has no zeros 
then this function corresponds to the ground state. 
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